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Abstract. We consider the evolution of a tight binding wave packet propagating in a time 
dependent potential. If the potential evolves according to a stationary Markov process, we 
show that the square amplitude of the wave packet converges, after diffusive rescaling, to a 
solution of a heat equation. 

1. Introduction 

It is generally expected that over long times the amplitude of a wave propagating in a 
weakly disordered background will be well described by a diffusion, at least in dimension 

> 3. This expectation stems from a picture of wave propagation as a multiple scattering 
process. Scattering off the disordered background results in random phases, and the build up 
of these phases over multiple scattering events leads eventually to a loss of coherence in the 
wave. Decoherent propagation of the wave may be understood as a classical superposition of 
reflections from random obstacles. As long as recurrence does not dominate the evolution, 
the central limit theorem suggests that the amplitude is given in the long run by a diffusion. 

So far it has not been possible to turn this heuristic argument into mathematical analysis 
without restricting the time scale over which the wave evolution is followed as in [21 El II]- 
One major obstacle is a lack of control over recurrence: the wave packet may return often 
to regions visited previously, denying us the independence needed to carry out the central 
limit argument. Indeed, the phenomenon of Anderson localization indicates that under 
appropriate conditions recurrence can dominate and cause complete localization of the wave 
packet. (It is worth noting that, since random walks are highly recurrent in dimensions d = 1 
and 2, the above heuristic analysis does not support diffusion in = 1 or 2, dimensions in 
which localization is proved {d = 1) and expected {d = 2) to dominate at any disorder 
strength.) 

A natural way to avoid recurrence difficulties is to bring a time dependence into the 
disordered background — we suppose that the environment evolves as the packet propagates. 
Here we consider a stochastic environment evolving independently of the wave packet. A 
natural assumption in this context is that the background changes in time according to a 
stationary Markov process. Such evolution equations have been proposed as an effective 
model for the propagation of wave packets in optical fibers [6J. 
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We consider here the simplest example of such a wave equation, namely the tight binding 
Markov- Schrodinger equation 



^0 e f{Z'^). 



where T is a translation invariant hopping operator on £^(Z'^), u(t) is a Markov process 
on a measure space Q and Wj, : f2 — M are measurable functions on Q. An elementary, but 
important, observation is that, so long as the time dependent generator H^(^t) = T+\vx{uj{t)) 
is uniformly bounded in time, the non-autonomous problem (11. ip has a unique solution ipt 
for any initial condition ipo G £^(Z'^), given for instance by the norm convergent series 

oo „ 

(1-2) = ^0 + '^i-'^y / H^iri) ■ ■ ■ H^{r„)^odri ■ ■ ■ dr„, 

n=i >^r„(t) 

with Tn(t) = {(ri,-- - ,r„) : < r„ < ■ ■ ■ < ri < t}. The evolution is easily seen to 
be unitary, WiptW = llV^oll- Thus a sufficient condition for solutions to (II. ip to exist is that 
||T|| < oo and sup^ sup^^ l^i:!'^)! < oo. 

We examine diffusion of the wave packet by considering the mean square amplitude 
E {\tpt{x)\'^) , where E (■) denotes averaging with respect to the random paths of the Markov 
process and also an initial distribution for uj{0). Diffusion is characterized by changes in 
position that scale as the square root of the elapsed time. Thus it is natural to look at the 
mean square amplitude in the scaling t i— > rt, a; i— > a/tx for a large parameter r. Since x 
is a discrete variable, to accomplish this rescaling we need to convolve E (|-?/'t(a;)p) with a 
function on W^. That is, we look at 

(1.3) A{x) = J2u{x-0^{\M0n 

with u a fixed "bump function" centered at 0. Let us suppose u> and / udx = 1, so that 
At{x) > and J At{x)dx = \\ipo\f- We interpret diffusion for the mean square amplitude 
as weak convergence of At{x) under diffusive scaling to a solution of a heat equation. That 
is, for suitable test functions 0, 

(1.4) / (f){x)T^ AT-t{\/Tx)dx '^~^°°> II^Aoll^ / — je~y^dx 



with D > 0. A sufficient condition for (11.41) . which requires no choice of a bump function, is 
obtained by a Fourier transform: 

(1.5) ^ ll^oire-*^!"^ VkGM'^. 



Following is a brief history of related studies. Ovchinnikov and Erikman obtained diffusion 
for a Gaussian Markov ("white noise") potential [7|. Fillet obtained results on transience of 
the wave in related models and derived a Feynman-Kac representation which we use below 
[8]. The evolution (11.11) was considered by Tchermentchansev P HD], who used Fillet's 
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Feynman-Kac formula to show that position moments such as l^l^-l^ (l"^*!^)!^) exhibit 
diffusive scahng, up to logarithmic corrections. More precisely, he obtained upper and lower 
bounds of the form 

(1.6) ^^(h^ < X^|xrE(|^i(x)|^) < tH\ntr\ t^oo. 

In the present paper we obtain diffusion ( II. 5p and show that (11. 6p holds for p = 2 with = 
z/_ = 0. While completing this manuscript, we learned of recent work of De Roeck, Frohlich 
and Pizzo on diffusion for a quantum particle weakly coupled to an array of independent 
free quantum fields [TT] . 

In the next section we state technical conditions which allow us to derive (II. 5p . These 
conditions are quite general and cover a large number of models of the form (II. ip . However, 
it may be useful to have at least one example in mind. So, to guide the reader, we close 
this introduction with a simple example of a potential for which we can derive diffusion. We 
call this the "flip process." The state space Q of the Markov potential is just { — 1, 1}^ , and 
Vx{^) = evaluation of the x^^ coordinate. So at any time t, the potential Vx{uj{t)) takes only 
the values ±1. Now suppose the process Lj{t) is obtained by allowing each coordinate Vx{uj) 
to flip sign at the times ti{x) < t2{x) < ■ of a Poisson process, with independent, identical 
Poisson processes at each site x. For this potential, our result implies diffusion (11.51) of the 
wave amplitude. 

2. Statement of the main result: diffusion of the amplitude 
2.1. Assumptions. We make the following assumptions: 

(1) {Existence of the Markov process and invariant measure): We are given a topological 
space Q, a Borel probability measure fi, and a collection {Pq : a G fi} of probability 
measures on the path space V = taken with the cx-algebra generated by Borel- 

cylinder sets, such that 

(a) {Paths are right continuous and start at a): For each a G f2, with Pq, probability 
one, every path uj{-) is right continuous and satisfies uj{0) = a. 

(b) {The Markov property holds): For any measurable ^ C P we have 

/ F^(^t){A)dF^{uj{-)) = F^{Si\A)) 
Jv 

where St is the backward shift on V, Stuj{-) = uj{- + t), so 
S^\A) = M-) : u;{. + t)eA}. 

(c) {Invariance of fi): For any measurable A <zfl, 

/ Fa{uj{t) G A)d/i(a) = /i(A) 
Jn 
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Let Ea(-) denote expectation with respect to Pq, 

(2.1) E„(F(cu(-))) = / Fico{-))dF^{u{-)), 
and similarly 

(2.2) E(-) = / E,(-)dMa) 



which is expectation with respect to the probability measure P(^) = f^Fa{A)dfj,{a) on V. 
By the invariance of /i under the Markov process we have 

(2.3) E(/(^(t)))= / /(a)d/x(a) 

Jn 

for any t and any / G 1/^(17). 

The Markov property, invariance of /i, and right continuity of paths show that 

(2.4) SJ{a) = E,(/(c.(t))) 

defines a strongly continuous contraction semi-group on L^(f2) (also on Lp{Q) for any 1 < 
P < oo). Indeed, the Markov property clearly shows this is a semi-group, and from the 
definition we have 

(2.5) {f,Stg)LHn) = E (/(cu(0))^?(a;(t))) , 
so by Cauchy Schwartz and 02.31) 

(2-6) \{f,Stg)LHn)\ < II/IIl^cc) II^IIl^ch) 

from which it follows that ||5't||^2^^2 < 1. The right continuity of the paths under now 
shows that St is strongly continuous, since any / G L'^{Q) may be approximated by bounded 
continuous functions and for bounded continuous f : Q —>■ M. we have 

(2.7) \\Stf-f\\l^n) = E{\f{u^{t))-f{umf) ^ 

by dominated convergence. 

The adjoint Sj of St is also a strongly continuous contraction semi-group, given formally 

by 

(2.8) Sjfia) = E(/(^(0))|^(t) = a), 

where the r.h.s. is a conditional expectation. Of particular importance to us is the generator 
B of sj defined by 

(2.9) Bi: = hm - (^j - Sj^p) 

t-*o+ t \ / 

on the domain T>{B) consisting of ip such that the limit on the r.h.s. converges in norm. 
The generator B is maximally dissipative, meaning Iie{ip, Bip) > for ip G T>{B) and no 
extension of B has this property. It follows that the spectrum of B is contained in the closed 
right half plane {z : Re^; > 0}. The adjoint i?^ of B is the generator for St- 
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Note that both St and Sj satisfy 

(2.10) Sd = Sjl = 1, 

where 1 denotes the function equal to one everywhere on Q. It follows that 1 G T>{B) and 
1 G V{B^) and that 

(2.11) Bl = Bh = 0. 

The orthogonal complement of 1 is the space of mean zero functions, 

(2.12) Llin) = |/GL2(fi) : ^ /(a)dM«) = o| . 

From fl2lTD it follows that L^(fi) is an invariant subspace for B and i?^. We require that B 
is strictly dissipative on this space: 

(2) {Gap condition and sectoriality of the generator): There is T > such that if / G 
V{B) and f{a)dfi{a) = 0, then 

Re{f,Bf)L^n)>^ [ |/(a)|' d/i(«). 

In addition, we require that B is sectorial, namely there is 7 < oo such that 
|lm(/,i?/)i.(f,)| < 7Re(/,S/)i2(^) 

for all / G V{B). 

One consequence of the sectorial condition on the generator is that a precise meaning can be 
given to the formal relation Sj = e~^^ using the Riesz functional calculus — see [H Chapter 
II, Section 4]. 

Finally, we require translation invariance for the hopping operator T, the Markov process, 
and the potential fa(a): 

(3) {Translation invariance of the hopping terms): T is a translation invariant hopping 
operator on £^(Z'^), 

TiIj{x) = ^h{x -y)i){y), 
y 

with h a function such that 

(a) {Self adjointness ofT) For every x, h{—x) = h{x)*. 

(b) {Non- degeneracy of T) For each non-zero vector k G M'^, there is some x G Z'' 
such that h{x) 7^ and k • x 7^ 0. 

(c) {Smoothness of the symbol) "^^Ixl"^ \h{x)\ < 00. 

It follows from (c) that h{'k) = e~''*^'^/i(x) is a function on the torus T'^. In 
particular, T is a bounded operator with ||7'||£2(2d)^^2(2d) = maxk |/i(k)|. 

(4) {Translation invariance of the process and invariant measure): There are //-measure 
preserving maps cTa, : f2 — > f2, x G Z"^, such that 
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and 

where % : V ^ V is the map %.{uj){-) = ax{uj{-)). 
(5) {Translation covariance, boundedness and non-degeneracy of the potential): The func- 
tions f : i7 ^ M are bounded, translation covariant 



Vx = vqo 

mean zero 

Vx{a)dfi{a) = 0, 



/ 



n 

and there is x > such that for all x,y E Z'^, x ^ y, 

(2.13) \\B-\v,-Vy)\\^,^^^>x. 

Since the Markov process is translation invariant, B commutes with the translations Txf{a) = 
f{axOi) of L^(fi). Thus ( I2.13P is equivalent to 

(2.14) \\B-\v.,-v,)\\^,^^^>x. 
for all X e Z'^, X ^ 0. 

The condition J f^d/i = can always be obtained by putting the mean of Vx into the 
diagonal part of the hopping term. Likewise, by absorbing the normalization into the disorder 
strength A, we may assume that ||'y2:||^2(f^) = 1. 

A very general class of models, which includes the flip model described above, is obtained 
by taking VL = S'^ for some set 5* C M and supposing that each coordinate a;(a;) of a; e f2 
undergoes an independent Markov process, with the processes at different sites identically 
distributed. We then set Vx{oj) = uj{x). In this case, the generator B is the sum of the 
individual generators for the processes at each site — more precisely the Friedrichs extension 
of that sum defined on the domain of functions depending on only finitely many coordinates. 
The above conditions are easily translated into conditions on the individual generator of 
the Markov process on S for each coordinate uj{x). For these models, the condition fl2.13l) 
is trivial since, by the independence of different coordinates and translation invariance, we 
have 

(2.15) \\B-\vx - Vy)\\l,^^^ = \\B-'vx\\l,^^^ + W^'Wllnn) = ^ ll^"'^o||'2(^) , x ^ y. 

2.2. Main result. The wave function ip satisfies a linear equation, but the square amplitude 
is quadratic in ip. To obtain a linear equation for the evolution of we consider the 
density matrix 

(2.16) pt{x,y) = Mx)Myr, 

which satisfies the evolution equation 
(2.17) 

dtpt{x, y) = -i ^ h{C) [pt{x -C,y)- pt{x, y + ()] - M^it)) - Vy{uj{t))) pt{x, y). 
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Note that \ipt{x)f = pt{x,x). 

More generally, we may consider the evolution equation fl2.17l) to be the basic dynamical 
problem, with an arbitrary initial condition pQ{x,y). The natural setting is for po to be a 
density matrix, namely 

(2.18) Po G VM := {p : Z'^ X Z'^ ^ C : p is the kernel of a non-negative definite, 

trace class operator on £^(Z'^)} . 

By virtue of the unitarity of the evolution fll.ip . the space T>Ai is preserved by the evolution 
(I2.17p . as is the trace trpo = J2x Poi^^-^)- 

Under the assumptions outhned in the previous subsection, we have the following 

Theorem 1. Any solution to (12.171) with initial condition po G VA4 satisfies 

(2.19) lim y e-'^-"E(p,t(x,x)) = [tr po]e-*^^^ ^''^(^^''-'^^ 



T—rOO ■ 

X 



with Dij = Dij{X) a positive definite matrix. Near A = 0, Dij{X) has an asymptotic 
expansion 

(2.20) D,,iX) = ^{Dl + OiX)). 



If furthermore |x|^po(a;, x) < oo, then 



2.21) lim ^ V|x|'E(pt(a;,x)) = [trpo]VA,i 



t^oo t 

X I 

Remark. As the proof will show, (I2.19P holds also for an initial condition po which is the 
kernel of a non-positive definite trace class operator. 

3. Augmented space analysis 

3.1. Augmented space and Fillet's Feynman-Kac formula. The starting point of our 
analysis is a "Feynman-Kac" formula due to Fillet which expresses E {pt{x, x)) as a matrix 
element of a contraction semigroup on an augmented Hilbert space, 

(3.1) n := L2(Z^ X X Q), 
where Z*^ x Z'^ x i7 is taken with the measure M 

(3.2) / f{x,y,u)dM{x,y,uj) = / f{x,y,u)dp{u). 

We think of a vector \1/ G 7i as a "random density matrix," at least if \l/(-, -,00) is the kernel 
of a non- negative definite trace class operator for p almost every u. We also think of Ti. as 
the tensor products f{Z'^) f{Z'^) ® L'^{n) or f{Z'^ x Z'^) ® L'^{n), using the notations 

(3.3) [ij ® ^ (E) f]{x,y,u) = ij{x)^{y)fiu), ij,^ei\Z''), f e L^Q), 
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and 

(3.4) [p®f]{x,y,uj) = p{x,y)f{u), p e f (Z'' x Z'') , f E L\Q). 
The Feynman-Kac-Pillet formula basic to our work is 

(3.5) npt{x,y)) = l,e-*(^^+^^^+^Vo® 
where 

(3.6) K^{x,y,u;) = J^HO ["^{x - C,y,uj) - ^{x,y + C,uj)] , 

( 

(3.7) V^ix^y^u) = {v^{uj) - Vy{uj))-^{x,y,uj), 

and the Markov generator B acts on 7i as a multiplication operator with respect to the first 
two coordinates, that is 

(3.8) B[p® f] = p®{Bf), pef{Z'^xZ'^),feL^{Q). 
In particular, we have 

(3.9) E{pt{x, x)) = {6^ O 5^ ® 1, e-'^po ® l)m 

where L = iK + iAl^ + B. This equation relates the mean square amplitude of the time 
dependent dynamics (12.171) to spectral properties of the non-self adjoint operator L. 

3.2. Fourier Transform. To perform a spectral analysis of L, it is useful to note that L 
commutes with a group of translations on 7i — a fact which encodes the distributional 
invariance of (11. ip under translations. Specifically, if we let denote a simultaneous shift 
of position and disorder, 

Si^{x, y, uj) = ^(x - ^, y - ^, a^uj), 

then we have 
Proposition 2. 

S^K = KS^, S^V = VS^, and S^B = BS^. 

Proof. The first two identities follow directly from the definitions of K and V; the last follows 
from the assumed translation invariance of the measure p. □ 

As K, V and B commute with a representation of the additive group Z^, we may simul- 
taneously partially diagonalize them by a Fourier transform. In the present context a useful 
transformation is the following unitary map: 

(3.10) $(x, uj,k) = J2 e-'^'-^S^^ix, 0,uj) = J2 e-^^-«^(x - ^, a^co), 

• : L2(Z2'^ xil)^ L^{Z'^ xVLxTd) 
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where = [0, 2nY is the d torus. One may easily compute that 

(3.11) K^{x, 00, k) := K^{x, cj, k) = ^ h{C) \^{x - C, oo, k) - e-'^<^{x - C, a^cj, k) 

C 

(3.12) V$(x, u, k) := V^{x, lu, k) = {v^^{uj) - voiuj))^{x, u, k) , 
and 

(3.13) m{x, u, k) = B^{x, uj, k), 

where B is understood to act as a multiphcation operator with respect to x G Z"' and k G T'^. 

The operators K, V and B act fiberwise over the torus T'^ — i.e., they act as multiphcation 
operators with respect to the coordinate k. Thus, eq. fl3.5p may be transformed into 



(3.14) E{pt{x,y))= / e-'^-^(4-,®l,e-*^''po;k®l)L2(z<*xQ)dW, 

where i denotes normalized Lebesgue measure on the torus T"^, 

(3.15) po;k(x) = ^e"*-Vo(x-y,-y), 

y 

and 

(3.16) Lk := iKk + iXV + B 
with V(f){x,uj) = {vx{uj) — vo{uj))(j){x,uj) and 

(3. 17) i?k0(x, u;) = Y, h{0 [^{x -Ccu)- e-^''<0(x - C, a^co)] . 

C 

In particular, 

(3.18) J2e-'''-^E{pt{x,x)) = {5o0l,e~''^^po.,k0l)mz''xny 

X 

This equation is the starting point for our proof of Theorem [H It indicates that the diffusive 
limit on the l.h.s. of fl2.19p can be studied via a spectral analysis of the semi-group e"*^'' for 
k in a neighborhood of 0. 

4. Spectral analysis of Lk and the proof of Theorem [T] 

In this section inner products and norms are taken in the space L'^{7/ x Vt) unless otherwise 
indicated. 
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4.1. Spectral analysis of Lq. To begin, let us consider Lk with k = 0. A preliminary 
observation is that 



(4.1) 



Lo^o ® 1 = 0. 



Ultimately, this identity is a consequence of the fact that E (tr pt) is constant in time. 

Let Pq denote orthogonal projection of L'^ilJ^ x Vt) onto the space Ti^ = ^^(Z'') {1} of 
"non-random" functions, 

(4.2) Po'^ix) = I ^(x,cj)d/i(cj). 

Jn 

Then Pf^ = (1 — Pq) is the projection onto mean zero functions 

(4.3) = i^^ix,u) : j ^ix,Lu)dfxiu) = oY 

The block decomposition of Lq with respect to the direct sum Hq © Hq is of the form 

(4.4) Lo = f ° ^ . 
^ ' \i\VPo iKo + B + iXP^VP^J 

Indeed, it follows from the definition of Kq that 

(4.5) KoPo = PoKo = 0, 
and we have seen in ^that 

(4.6) PoB = BPo = 0. 

(The identity PqB = follows since -B^Pq = 0.) Thus, ranPg-^ = Ti^ is an invariant subspace 
for iKo + B and iKo + B vanishes on ran Pq = Tio- Since the potentials Vx{uj) are mean zero, 
PqVPq = and gll) follows. 

Using the block decomposition (14. 4p we now prove the following 

Lemma 3. For each A > there is 6\ > 0, with 6x > cA^ for A small, such that 

(4.7) a(Lo) = {0}US+ 
where 

(1) is a non- degenerate eigenvalue, and 

(2) ^+C{z : Re^ > 6x} ■ 

Proof. First note that ReLg = KeB > in the sense of quadratic forms. Thus by the 
sectoriality of B 



(4.8) lm($,Lo$) 



< 



Ko + \V 



+ |lm($,P$)| < ||/i||oo + 2A + 7Re($,Lo$), 
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if ||<l>|| = 1. It follows that the numerical range of Lq, Num(Lo) = |($,Lo$) : ||$|| = 1 
is contained in 



(4.9) 



A/+ : = < z : Kez > and I Im 2; I < 



+ 2A + 



7 Re 2;| . 



Since (j(Lo) C Num(Lo), we may restrict our attention to 2; G 
Now fix 2; G and consider the equation 



(4.10) 



(^0 - z) 



® 1 \ ^ / -z iXPoV \ f(j)®l 

$ / \\\VPn Pj-Lr.Pj- - 7} \ $ 



® 1 



V^v " ^ -^O-i ^/ 

for (0 ® 1, $) G 'Ho © ^0 given {i/j ® 1, ^') G © Hq • By the gap condition on 5, 



(4.11) 



so the second of the two equations (14.101) may be solved provided Re ^ < ^ to give 

(4.12) $ = {P^LqP^ - zY^Ul -iXVcj)®! 
Thus the first equation of (14.101) reduces to 

(4.13) [T{z)- z](P(^l = PoV{P^LoP^ - z)~H, 
with 

(4.14) r(^) = \^PoV{P^L^P^-z)-^VP^. 

Thus Lo — z is boundedly invertible, for Rez < ^, if and only if z ^ a{T{z)). However, 

(4.15) Re(0® l,r(z)0® 1) 

= A2(0 ® 1, V{P^lIp^ - z*y' (Re 5 - Re;z) {PqLoP^ - z)-^V(f) ® 1) 

2 



T 



Re 2 



-1 



A^ 



[^-Re-}j {b-^P^{Lo-z)P^Y^ B-^V<P®1 



where the inverse of B is well defined since Vcj) ® 1 G Ti^ = ranPj^. Furthermore B ^ 
bounded on H^, with ||P"^Po^|| < T. Thus B-^P^{Lo - -2)^0^ is bounded, 



(4.16) 



B-'P,'-{Lo-z)P,^ 



<1 + T 



+ 2A + l^l 
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and it follows that 
(4.17) 

Re(0 ® 1, T{z)^ ® l)£2(zd) > 



Re 2; 



1+T \\h 



2A + 1^1 
1 



5-V0(g) 1 



1 + T i\\h\\^ + 2X + \z\ 



where we have made use of the non- degeneracy assumption on the potential. 
For z e A/V n (Re 2 < |; } we have 



(4.18) 
Thus 

(4.19) 



\h\ 



2\+\z\ < 2||/i||oo +4A + (1 + 7)-. 



Rerr^) > 



1 - T Re 2 



T 



Ax 



4rA 



:i - no) 



_2 + 7 + 2T||/i||, 

with IIo the projection of Hq onto Sq. It follows that 

(1) r(0) is invertible off the range of IIo, so is a non-degenerate eigenvalue of Lq, 

(2) r{z) - z is invertible if 2; G A/V, 2; 7^ 0, 

-\ 2 

Ax 



(4.20) 



r 1 — r 

Re 2; < — , and Re^; < — — — 
- T ' - T 



2 + 7 + 2T||/i|U + 4TA 



for some r G (0, 1). 
Optimizing the choice of r gives the following explicit expression for S\: 

1 A^X' 



(4.21) 



^ ^2 + 7 + 2T\\h\\^ + ATXV + AV 



□ 



The spectral gap 5x has consequences for the dynamics of the semi-group: 

Lemma 4. Let Qq = orthogonal projection onto 5o ® 1 x ^) • Then e^*^" ( 1 — Qo) is 

a contraction semi-group on ran(l — Qq), and for all sufficiently small e > there is > 
such that 



(4.22) 



,-tLoi 



1-Qo 



|ran(l-(3o) 



Proof. That e^*-^''(l — Qq) is a semi-group on ran(l — Qq) with generator Lq^ = ZqI 
is clear. Since cr{Lo) C {Re;z > 5a}, it is known that a sufficient condition for fl4.22p to hold 
is for e~*''^°(l — Qq) to be an analytic semi-group (see [H Chapter IV Corollary 3.12]). A 
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necessary and sufficient condition for e *'^°{l — Qo) to be analytic, given that Lq is maximally 
dissipative, is that the following estimate 



(4.23) 



< 



C 



ran(l-Qo) 



ImCI 



holds for all C G C with ReC < 0. (See [H Chapter II Theorem 4.6].) 

For any invertible closed operator A, we have IIA^"*^!! < 1/inf {|2;| : z G Num(y4)}, with 
Num(y4) the numerical range of A. Since Num(L[,°^) C 7V+, with ^^+ as in (14.91) . we have 



(4.24) 



< 



1 



ran(l-Qo) 



dist(c,A/;)' 



It follows that (14.23P holds for all ( outside a rectangle of the form R = [—a, 0] + i[—M, M]. 
On the other hand, by Lemma [3], this rectangle R is contained in the resolvent set of Lq\ 
Since R is compact, we have 



(4.25) 



sup 

Cei? 



which is stronger than (14.231) for ( & R. 



< oo, 



ran(l-Qo) 



□ 



4.2. Analytic perturbation theory for Lk. Now that we have established a strict spectral 
gap for Lq, it follows that the gap persists in the spectrum of Lk for k sufficiently small. 
Indeed, 



(4.26) 

so 

(4.27) 



[Lk - Lo\(f){x,uj) 



c 



■^)0(a;-C,cTc^), 



L\i — Lq 



< c\k\. 



Thus, an immediately corollary of Lemma [3] is 



Lemma 5. // |k| is sufficiently small, the spectrum of L]^^ consists of: 

(1) A non- degenerate eigenvalue E{k) contained in Hq = {z : \z\ < c|k|}. 

(2) The rest of the spectrum is contained in the half plane Hi = {z : Re z > 6\ — c\k\} 
such that ifo n iJi = 0. 

Furthermore, -E(k) is in a neighborhood ofO, 



(4.28) 



E{0) = 0, VE(0) = 0, 
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and 

(4.29) d,djE{0) = 2Re{diKoSo®l,[Lo]~^djKo6o^l) 

= 2Re J2 Xiyjh{^h{y){S,®l,[T{0)]-^Sy®l), 



with 

(4.30) r(0) = X^PoV{Po^LoP,'-)-'VPo 

as in the proof of Lemma\^ In particular, didjE{0) is positive definite. 

Proof. These are all standard facts from analytic perturbation theory — see for instance 
[5]. A sketch of the proof is as follows. First, it is a general fact that the spectrum moves 
no further than the norm of the perturbation, so 

(4.31) a(Lk) C {2 : \z\ < c\k\} U {z : Re z > 5x - c\k\} . 

For sufficiently small k the two sets on the r.h.s. are disjoint and we may fit a contour C 
which winds around the origin between them. The (non-Hermitian) Riesz projection 

(4.32) = J- f -Ji^dz 

27n Jc z-Li, 

is rank one at k = and continuous as a function of k. It follows that Qk is rank one as 
long as o"(Lk) does not intersect the contour C. Thus for small k the only spectrum of Lk 
in a neighborhood of zero is a non-degenerate eigenvalue, with associated eigenvector in the 
one-dimensional range of Qk- 

Let us call the eigenvalue -E(k), and the associated normalized eigenvector $k- Clearly 
E{Q) = and $0 = 5o ® 1- Since 

(4.33) LkQk = ^(k)gk 

we may compute the derivatives of E(k) by differentiating Lk and the projection. In partic- 
ular, V-E(k) is given by the so-called Feynman-Hellman formula 

(4.34) VE{k) = ($k,VLk$k), 

from which it follows that V-E'(O) = since VLk = iV-ft'k is off-diagonal in the position basis 
on 7^0- Similarly, we have 

(4.35) d,d,E{k) = ($k,5,a,£k$k) + (5^^k$k,(l-Qk)2k'(l-Qk)a,£k$k) 

+ (S.Lk'I'k, (1 - Q^)L^\l - Qk)9.Lk$k) 

The first term on the r.h.s. vanishes at k = and the other two combine to give the identities 
claimed in the Lemma. Since [r(0)]~^ is positive definite, it follows from the non- degeneracy 
condition (3b) that didjE{0) is positive definite as well. □ 

Again, we obtain dynamical information about the semi-group e"*^'': 
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Lemma 6. If e is sufficiently small, then there is < oo such that 



(4.36) 

for all sufficiently small k. 



Proof. Since Num(Lk) C A/^, with A/+ as in (14. 9 p in the proof of Lemma[3l this is essentially 
identical to the proof of Lemma HI For k in a compact neighborhood of the origin, we can 
choose the bound Cf uniform in k. □ 



4.3. Proof of Theorem [T]. A first observation is that it suffices to prove diffusion (12.191) 
under the assumption that the initial density matrix po satisfies 

(4.37) ^|po(a;,i/)| < oo. 

To see this, it is useful to note: 

(1) The evolution (I2.17P preserves the trace norm of Pq. 

(2) Any p G VJ^ may be approximated, to arbitrary precision, by an operator po G VJH 
satisfying (14.370 . 

In more detail, recall that the trace norm of an operator A is 

(4.38) Pll^^ = sup{|trA5| : 5 is finite rank and < 1}. 

Since the evolution (I2.17P is given by pt = UpoW, with U = U(t,0) the unitary propagator 
of (II. ip . we see that ||pt||^^ = ||po|lrr ^n the other hand, given po of trace class and e > 
we can find po satisfying (I4.37P and such that ||po — PoIIti < ^- Indeed, we may take 

po{x,y) if|x|,|?/|<L 
if \x\ > L or \y\ > L 

with L sufficiently large. Then, 



(4.39) 



(4.40) 



J^e '^■''E{pt{x,x) - pt{x,x)) 



< J2^i\ptix,x) - pti^,x)\) < E (Ilpi - MrJ < 

X 

If diffusion (12.190 holds for any po satisfying (I4.37P we learn that 



(4.41) 



lim sup 



^e 'v^'''^E(p^t(x,x)) - (trpo)e ^ 



< 2e. 



which gives diffusion for po in the limit e ^ 0. 

Turning now to po which satisfies (I4.37P we see that 

(4.42) Po;k(x) = 5^po(x + y,i/)e''^-2' 
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defines a function tliat is uniformly bounded in P'ilJ^^ as k varies tlirougli tlie torus 



(4.43) 



< 



XllPo;k(a;)| < ^|po(a;,l/)| := M < oo. 



x,y 



By (Km we liave 

(4.44) Ve-'7?^-"E(p.i(x,a;)) = {60 ^ l^e-^'^'^/^p,.^^® 1). 



Letting Qk denote tlie Riesz projection onto tlie eigenvector of Lk near zero — see (14.321) 
in the proof of Lemma [5] — , we have 

(4.45) Ve-'V?i^-"E(p.i(x,x)) = e-^*^(^/^^)(5o ® 1, Q^kPo-^k ® 1) 

X 

+ (5o®l,e-^*^''/v^(l-Q 1 k)Po; 1 k®l), 

for r sufficiently large. By Lemma [6l the second term in f l4.45p is exponentially small in the 
large r limit, 



(4.46) {So ® 1, (1 - Q 1 k)e-"*'^''/-^Po; 1 k ® 1) 



< 



(1 - Qi^)e-^'''-/^ 



PO;^k ® 1 



< MCe"^*^'^^"'"''''!/^) ^ 0. 



Regarding the first term in ( 14.451) . we have by Taylor's formula, 

(4.47) ^(k/v^) = ^5^9,a,E(0)k,k,. + o(^), 

since E(0) = V^(0) = 0. Thus 

(4.48) g-rtE(k/v^) ^ g-4E.., 9,9,i?(0)k,:k, ^^(^^^)_ 

Putting together (OHj) and (06|) yields 

(4.49) ye"'^'^-"E(p,i(x,a;)) = e"*^ ^'^^^^'^''^^ ((5o ® 1, p^.^k ® 1) + o(l) 



= e-*5E-^-^^^(°)'''''^Po.^k(0)+o(l) ^ e-*5A..,«.«.^Wk,k,^^^^(Q) 

since Qj^^o ®1^5o®lask^O and Po;k(0) is continuous as a function of k. 

Eq. (I4.49P completes the proof of diffusion (12.191) . with the diffusion matrix given by 
Di j = ^didjE{0). From the explicit expression (14.291) for didjE{0) in Lemma [5], we see that 
the asymptotic form ( 12.201) holds, with 



(4.50) = Re (diKoSo ® 1 



-1 
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To derive (12.211) . diffusive scaling for ^ ^))' ^ote tliat by (13.181) 

(4.51) 



k=0 



Expanding tfie r.li.s. we obtain 

(4.52) jJ2x^E{pt{x,x)) = -y(5o® l,e-*^°AkPo;k|k=o® 1) 



-- y"(5o ® 1, e-'^^ di Po;klk-o ® 1) 



i=l 



II 



k=0 



P0;0 



III 



Tlie first term is negligible as t ^ oo, 
(4.53) 



I = -\{So®l,e *^°AkPo;k|k=o® 1) = -jJ2^^Po(^'^^ = ^(7 



t 



since e ^^°6q ® 1 = 5o ® 1- 

To proceed, we need to recall the formula for differentiating a semi-group, 



(4.54) 



Thus the second term on the r.h.s. of (14.521) is equal to 



(4.55) 



II 



t 



I V(5,Lo5o ® 1, e-^^°(l - Qo)dipo.,o ® l)ds, 
i=i 



where we have recalled that 



(1) 



t 



(2) e~'^"Qo = Qo, and 

(3) Qod^Lo = 0. 

By Lemma m we see that II is negligible in the large t limit: 



(4.56) 



nil < 1 re-'^'^-^Ms = o(-]. 
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It remains to compute the third term III, for which we have 

(4.57) III = -- / (AkLk|k=o5o®l,e-^^°(l-Qo)po;o®l)ds 
t Jo 



Ilia 



+ 7 r rZl(Wo®l,e-(^-^)^ni-Qo)a«2oe-^^°(l-Qo)po;0®l)drds 
t Jo Jo 



i=l 

V — 

Illb 



Iff V(5,Lo5o ® l,e-(^-'-)2''(l - Qo)5.^ogo Po;o ® l)drds 



IIIc 

Similar to what we had for II, we find that Ilia = O (j), and also that 

(4.58) |IIIb| <jj^j^ e-^(^^-^)drds = O . 
Finally, for IIIc we have 

(4.59) IIIc = ^(9,Lo5o®l,2o'(l-Qo)5i2ogoPo;o®l) + oQ^, 

^ ^ ^ /I 

= (trpo) Y^{d,Ko6^ ® 1, L-'d,Ko6o ® 1) + O f - 

since VL^ = iV-ft^k- Comparing with the expression fl4.29p for didjE{0) yields (12.211) . □ 
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